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ABSTRACT 

We present a quantitative analysis of the largest contiguous maps of projected mass 
density obtained from gravitational lensing shear. We use data from the 154 deg^ cov- 
ered by the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Our study 
is the first attempt to quantitatively characterize the scientific value of lensing maps, 
which could serve in the future as a complementary approach to the study of the dark 
universe with gravitational lensing. We show that mass maps contain unique cosmolog- 
ical information beyond that of traditional two-points statistical analysis techniques. 

Using a series of numerical simulations, we first show how, reproducing the 
CFHTLenS observing conditions, gravitational lensing inversion provides a reliable 
estimate of the projected matter distribution of large scale structure. We validate our 
analysis by quantifying the robustness of the maps with various statistical estima- 
tors. The same process is then applied to the CFHTLenS data. It is found that the 
statistical properties of the projected mass are fully consistent with the cosmological 
analysis performed on the shear signal discussed in CFHTLenS companion papers. 
The maps also lead to a significant measurement of the third order moment of the 
projected mass in agreement with analytic predictions, and to a marginal detection 
of the fourth order moment. Tests for residual systematics are found to be consistent 
with zero for the statistical estimators we used. A new approach for the comparison 
of the reconstructed mass map to that predicted from the galaxy distribution reveals 
the existence of giant voids in the dark matter maps as large as 3 degrees on the sky. 
Our analysis shows that lensing mass maps not only is consistent with the results 
obtained by the traditional shear approach, but also enable new analysis techniques 
such as peak statistics and the morphological analysis of the projected dark matter 
distribution. 

Key vifords: dark matter - large-scale structure of the Universe - gravitational lensing 



2 L. Van Waerbeke et al. 



1 INTRODUCTION 

Gravitational lensing is a powerful tool for the study of 
the dark matter distribution in the Universe. The statis- 
tical analysis of the distortion and magnification of lensed 
galaxies contains unique information on structure formation 
processes. For example, two point statistics of the shear or 
convergence can be used to constrain the dark matter power 
spectrum, the growth of structure, and cosmological param- 
eters (see |Munshi et~aL||2008| [Hoekstra fc Jain, 2008, , for 
recent reviews). Over the past ten years, multiple groups 
have reported improved constraints on the mass density pa- 
rameter and the power spectrum normalisation using the 
shear correlation function. However, it is well known that 
gravitational lensing contains a lot more information than 
the amplitude and shape of the mass power spectrum. In 
particular the distortion (shear) can be used to reconstruct 
the projected mass density to make mass maps. 

The purpose of this paper is to explore to what extent 
mass maps can access the additional information not cap- 
tured by two points statistics. Several groups have explored 
different theoretical routes beyond two points statistics. For 
instance, higher order shear measurements are a sensitive 
measure of the gravitational collapse process through mode 



coupling in the non-linear regime (Bernardeau et al. 1997; 
Van Waerbeke et al.|200Tl|Takada fcJ ain 2003; Kilbingcr & 
Schneider|2005 1, and can also be used as an indicator of non- 
Gaussian non-linearity in the primordial dark matter distri- 
bution ([Takada fc Jain|[2004l [Valageas et al.||2005[|. More 



exotic statistical estimators involve global statistical tools, 
such as the Minkowski functional, as opposed to local mea- 
surements based on the shear correlation functions. Mor- 
phology of large-scale structures for instance ( [Mecke et al.| 
1994 1, is also a probe of the non-linear processes in action 
during structure formation (Sat o et al.|2001 |. 

The statistics of peaks is also an estimator that can be 
used to probe how many dark matter haloes there are at 
some redshift within some mass range (|Jain fc Van Waer^ 



beke|2000||Kratochvil et al.|2010||Maturi et al.|2010[), and it 



can potentially constrain the halo mass function ( Yang et al 
|2011[ ). The effect of the large-scale dark matter environment 
on the galaxy-galaxy lensing signal could be investigated; for 
instance if some well known galaxy property (e.g. star for- 
mation rate, dust content, stellar population) is dependent 
on the average dark matter content well beyond the galactic 
halo scale, mass maps can be used to perform this measure- 
ment as function of the projected density. 

Most of these non-local statistics cannot be expressed 
as a combination of a subset of moments of the shear, which 
are, by definition, local measurements. Global, or non-local, 
statistics (such as peaks and morphology) require access to 
the full 2-dimensional mass maps obtained from the non- 
local gravitational lensing reconstruction. We note that dark 
matter maps can also be used to understand global features 
in the baryon/dark matter relation via cross-correlation with 
maps from different wavelengths surveys, whether Sunyaev- 
Zel'dovich, X-ray, or atomic hydrogen. The use of mass maps 
is therefore complementary to the analysis of the shear and 
magnification two-point correlation functions. This moti- 
vates the need to have reliable mass maps from gravitational 
lensing data. 

Mass reconstruction is a common tool for the study of 



galaxy clusters ( Hoekstra et al.||2000 Clowe et al.||2006a|b 



Heymans et al.|2008 1; it is mostly used to qualitatively eval 
uate how well the baryonic and dark matter distributions 
match, and to explore the possible existence of dark clumps. 



A few rare cases suggest their possible existence ( Erben et al 



2000| [Mahdavi et al.||2007| |Jee et al.||2012[ ), although the 
significance of these detections is currently unclear. Mass 
reconstruction in cluster environments concern regions of 
relatively high lensing signal and small angular scale. Mass 
re construction beyond galaxy cluster scales was performed 



in |Massey et al.| ( |2007[ ) with the 1.64 deg^ area of the COS- 
MOS survey. However, our work shows that sampling vari- 
ance is still important at even larger scales, implying that 
a relatively small survey like COSMOS does not provide a 
fair sample of the dark matter distribution. 

In this paper we use the 154 deg^ CFHTLenS lensing 



data jErben et al.|20l"2|jHeymans et al.|2012b[ [Miller et al 



2013 [Hildebrandt et al.|2012 1 to perform dark matter mass 
reconstruction over angular scales of several degrees. We do 
not attempt to perform a 3-dimensional mass reconstruction 



(,Massey et al. 2007 


Simon et al. 2012 VanderPlas et al. 


2011| [Leonard et al. 


|2012| because such reconstruction us- 



ing ground-based data would only yield meaningful results 
for galaxy clusters relatively close in redshift, z < 0.3, or for 
very massive and rare clusters at higher redshift, z < 0.6 
(1 Simon et al. 20091. We explore the measurement of con- 



vergence statistics from gravitational lensing maps and the 
connection between dark and baryonic matter up to a few 
degrees in angular scale. In particular we focus on the im- 
portance of the effect of noise in the reconstruction, and how 
to account for its effect on the measurement. To this end, 
ray-tracing simulations are used to test our reconstruction 
and analysis procedure. 

This paper is organized as follows. Notation and mass 
reconstruction strategies are described in Section 2. The 
analysis of reconstructed mass maps based on N-body sim- 
ulations is discussed in Section 3. In Section 4 we show the 
results from the CFHTLenS data. Section 5 presents our 
new approach for the comparison of dark matter and bary- 
onic matter, with an emphasis on the detection of voids. We 
conclude in Section 6. 



2 MAP MAKING AND COSMOLOGY 
2.1 Mass reconstruction 

Mass maps are proportional to the projected mass convolved 
with the lensing kernel. They can be constructed from the 
shear measurement, the details of this process and the no- 
tations are given below. 

A galaxy at position 6 = {61,02) on the sky, is char- 
acterised by its redshift z and its shear components "/i{0), 
where i = 1, 2. In the following, we always work on the flat 
sky tangent plan. The relation between "f and the conver- 
gence K involves the gravitational deflection potential '5(0) 
defined as: 



7,(0) = G,*(0), 



(1) 



where Gi = 9ii — 922 and G2 = 29i2. The function is 
given by a line-of-sight integral of the 3-dimensional matter 
gravitational potential $ via: 
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^ ' Jo fK{w')fK{w: 



^{fK{w')e,w') 



(2) 



The comoving radial distance at redshift z is given by w{z), 
and fxiw) is the corresponding angular diameter distance. 
The comoving distance at the source redshift is Ws- The 
3-dimensional gravitational potential depends on the mass 
density contrast 5 via the Poisson equation: 



3-ffo ^0 . 
2a ' 



(3) 



where Ho is the Hubble constant, f2o is the matter density, 
and a is the scale factor. The convergence k,{6) is defined as 
the line-of-sight projection of the 2-dimensional Laplacian 
of t/j: 



= A*(0), 



(4) 



where A = flu + 922, and the derivative 833 of the gravita- 
tional potential along the line-of-sight vanishes on average 
due to the Limber approximation (Kaiser 19981. The ob- 
servable quantity is the reduced shear g: 



g = 



7 

1 - k' 



(5) 



Due to the intrinsic ellipticity of galaxies, galaxy shapes pro- 
vide a noisy estimate of This estimator, noted e°''°, is 
given by : 



1 - ge'" 



(6) 



The intrinsic ellipticity e'"' is a generic term which contains 
a contribution from the intrinsic galaxy shape and from the 
measurement noise. The ensemble average (e°'"') is an unbi- 



ased estimator of the reduced shear g ( Kayser & Schramm 



19881. In the case of CFHTLenS, a complete description of 



the measurement noise is given in Heymans et al. (2012b I 



and |Miller et al.| ( |2013) . 

Weak lensing studies assume that both the shear 
and convergence amplitudes are much smaller than unity, 
^ 1) so that one can use the weak lensing ap- 
proximation (g ~ -y) for mass reconstruction. In that case 
Eqje] reduces to e"''" ~ 7 -I- e'"* , and it follows that by 
comparing Eq[T] and Eq|4] one can reconstruct the conver- 
gence from the observed shear in Fourier space (Kaiser & 
Squires|1993[ ). For the weak lensing approximation, we first 
write the relation between the shear and convergence Fourier 
components k(s) and 'y(s), where we define the wavevector 
s = (si,S2) = 2ii/6 as the 2-dimensional analogue of 6 in 
Fourier space: 



k(s) = 



^1 ~ ^2 



7i(s) + 



p , p 72(s). 

n/-j^ ~I~ '■'2 



(7) 



Following Kaiser & Squires (19931, the shear data are first 



regularised with a smoothing window and then EqjT] evalu- 
ated from the smoothed Fourier components. Since real data 
also contain masked regions (due to bright stars, and other 
aberrations in the image), the smoothing takes into account 
the number of pixels being masked within each smoothing 
window, so that the resulting averaged quantity in that win- 
dow is not biased. The practical procedure for obtaining a 
mass map from reduced shear is the following: 



• The data e°''=(0i 



e° are first placed on a regu- 



lar grid 6ij and then smoothed. At pixel location 8ij, the 
smoothed ellipticity e^j is given by: 

where Woq (6) is a normalised Gaussian smoothing window: 



exp 



\9l 
291 



(9) 



and w(6) is the weight associated to the measurement noise 
of e{6). Details on the weight specific to CFHTLenS data 
are given in Section 4.2. 

• The reconstructed mass is obtained from Eq[7] Note 
that the Gaussian filter has infinite spatial extension, there- 
fore none of the pixels in the final grid have a value of zero. 
There are however pixels with larger noise (fewer objects) 
than others, especially when a mask overlaps with the cen- 
tral region of the filter. The noisiest pixels in the recon- 
structed map are removed if the effective filling /actor within 
the Gaussian window is below 50 percent. This aspect is 
discussed later in Section 3.2 where the noise in the recon- 
structed maps play an important role in the cosmic statistics 
measurement from the convergence. 



The alternative to the Kaiser & Squires ( 1993 1 mass re- 



construction approach is the full non-linear mass reconstruc- 
tion as described in Bartelmann et al. (19961. A detailed 
comparison between the two methods is left for a future 
study. For the purposes of this paper, the KS93 reconstruc- 
tion is sufficient as we are not trying to recover the lens- 
ing signal in the non-linear regime at sub-arcminute scale. 
We note however that the ability of the non-linear mass re- 
construction to recover large scale structure statistics was 



demonstrated in Van Waerbeke et al 
As shown in 



(19991. 



Crittenden et al. (20021, the shear can 



split into E and B modes, which correspond to the curl-free 
and curl shear components respectively. In the absence of 
residual systematics the scalar nature of the gravitational 
potential leads to a vanishing B-mode. As shown in |SchnePl 
(19981, the split between E and B modes can be 



der et al. 



made by applying the transformation (71,72) to (—72,71), 
which is the same as rotating each galaxy by 45 degrees. In 
the following, mass maps are reconstructed using the E and 
B modes and the corresponding convergence is called 
and Kb- As we will see, it will be necessary to distinguish 
between the reconstructed (noisy) map, and the underlying 
true convergence value. We therefore introduce three ad- 
ditional convergence terms. We call Kobs the reconstructed 
convergence such that: 



^obs — I^E ~t- ^r; 



(10) 



where k_b is the true underlying signal and Kran is the recon- 
struction noise. This equation is the convergence equivalent 
to Eq.(|6| in the weak lensing regime. We also introduce the 
convergence kx, reconstructed from the 45 degrees rotated 
galaxies such that: 



(11) 



where hb should be consistent with zero if the residual sys- 
tematics are negligible. The validity of this statement will 
be verified statistically. 
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2.2 Cosmology 

The noise-free convergence map Keg{6) smoothed with a 
Gaussian filter of size 9o can be expressed in terms of the 
3-dimensional mass density contrast S{fK{w){6),w): 



K-eo {0) 



2 



AO' 



1 



/o «(™) 

5{fK{w){e-e'),w)We,{e' 



(12) 



where Wsq (0) is the normalised Gaussian window given by 
Eq[9] and now wh represent the radial distance at infinite 
redshift. The function g{w) accounts for the distribution of 
source redshifts: 



g{w) = 



Aw ps{w ) — 

.fK(w') 



(13) 



where ps{w{z)) is the source redshift distribution. Eq.(12l 
represents the convergence map obtained from the recon- 
struction process after the shear data have been pixelated 
and smoothed as described in Section 2. Therefore the po- 
sition on the sky are pixels on the regular grid where the 
map has been computed. The second order moment of the 
convergence {K^{do)) can be expressed from the convergence 
power spectrum : 

{^')Bo = l^Y.^-2l + l)CtWn0o), (14) 

where Wi{Oo) are the multipole moments of the smooth- 
ing window. For the Gaussian window given by Eq.([9|, the 
multipole moments are given by Wi{6o) — exp{—l 0^/4). 
The power spectrum Cf can be derived from Eq[T2] assum- 
ing that the small angle and Limber approximation applies 
| |Miralda-Escudc 1991 Kaiser|1992[ ): 



cr = 



d»4Mp3, 



I 



fKiw)' 



(15) 



(16) 



fK{Ws) 

The 3-dimensional mass power spectrum is defined as: 

{5*(k,u;)5(k',«;)) = (27r)^fc(k-k')P3i5 ik-w). 

The third order moment {i<^)eo can also be expressed ana- 
lytically using the convergence angular bi-spectrum Bi^i^i^, 
the smoothing filter multipoles WiiOo) and the Wigner 3-j 
symbols: 



{2h + l){2l2 + l){2h + 1) 
4tt 



-y 

li h h 


X B^^i^t.^w,,{eo)Wi,{eo)Wi,{eo 

with the angular bi-spectrum given by 



Stt 



Sm't (2h + l)(2l2 + l){2l3 + l) f h h h 



An 







Aw 



d'^{w)fK{w) 



X BzD 



fK{w) ' fKiw) ' fK(w) ' ™ 



where Bsd is the three dimensional bi-spectrum. Ban is cal- 
culated in the non-linear regime using non-linear extensions 
of perturbation theory. Expressions for gravitational lensing 



are derived in Van Waerbeke et al. ( 2001 1, based on the work 



of Scoccimarro & Couchman (2001 1 



3 SIMULATIONS 

3.1 The CFHTLenS data 

The purpose of the following Section is to validate our mass 
reconstruction approach using mock catalogues that repli- 
cate the true CFHTLenS observing conditions. CFHTLenS 
spans a total survey area of 154 deg^, covered with a mo- 
saic of 171 individual pointings observed by the one square 
degree imager at the Canada France Hawaii Telescope. The 
survey consist in four compact regions called Wl, W2, W3 
and W4, which cover approximately 72, 36, 50 and 25 deg'^ 
respectively. Detai ls on the data reduction are described in 



Erben et al. 



( 2012 1. The effective area is reduced to 120 deg 



by the masking of bright stars, artificial and natural moving 
objects, and faulty CCD rows. The observations in the five 
bands u*g'r'i'z' of the survey allows for the precise mea- 
surement of photometric redshifts (Hild ebrandt et al.|2012 |. 
The shape measurement with lens&t is described in detail in 



Miller et al. (20131 



3.2 Mock catalogues 

The mock catalogues are constructed from a mixture of 
real data (the noise, galaxy position and masking structure 
of CFHTLenS) and simulated data (the shear signal from 
N-body simulations). The N-body simulations used in this 
work are described in Harnois-Deraps et al. ( 2012 I. The pro- 



cedure for generating mock catalogues is the following: 

(i) Projected shear and convergence maps are constructed 
from the combination of redshift slices for all source red- 
shifts available in the simulation The maps are sampled 
on a 1024 x 1024 grid and cover 12.4 deg^ each, which cor- 
responds to a native pixel size of 0.21 arcmin. These maps 
represent individual 'tiles'. The CFHTLenS fields Wl, W2, 
W3 and W4 are covered with the maximum number of non- 
overlapping tiles that fit within their respective areas. Given 
the are of the four CFHTLenS mosaic fields, one can only 
cover Wl and W3 with 4 tiles each, and W2 and W4 with 
one tile each. The final simulated area is therefore 124 deg^; 
this is similar the 154 deg^ covered by CFHTLenS. Note 
that each of the tiles used in our simulations are indepen- 
dent lines-of-sight, therefore we expect sampling variance to 
affect the mock catalogues less than the four CFHTLenS 
fields, where only four fields are statistically independent. 

(ii) The 2d positions on the sky of all galaxies in 
CFHTLenS are preserved. The orientation of each galaxy's 
ellipticity is randomized, but its amplitude is held fixed as 



1 The source redshifts are 0.025, 0.075, 0.126, 0.178, 0.232, 0.287, 0.344, 
0.402, 0.463, 0.526, 0.591, 0.659, 0.73, 0.804, 0.881, 0.961, 1.071, 1.215, 
1.371, 1.542, 1.728, 1.933, 2.159, 2.411, 2.691, 3.004 
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a.35 



Figure 1 . Mass reconstruction for one simulated line-of-sight cov- 
ering 12 square degrees. The continuous background map with 
white contours represents the reconstructed lensing mass (con- 
vergence) with masks shown by the black regions. The white con- 
tours show the 1, 2, 3 and 4 a contours, while the black contours 
show the 1, 2, 3 and 4 a levels in the noise-free map. a is the 
convergence rms measured on the mass reconstruction. 



this is used as the only source of shape noise in the mock 
catalogue. The redshift of each galaxy is resampled using 
redshift probability distribution obtained from the photo- 
metric redshift distribution function. We refer the reader 
to 



Benjamin et al. (20121 for the details. The galaxies are 



then placed on the mosaic tiles and each galaxy's redshift 
determines uniquely the simulated shear and convergence. 
This is calculated following Eqijl] and |4] with the reduced 
shear computed from Eq[5] The final 'observed' ellipticity 
in the mock catalogue is obtained by Eq|6] The masks are 
also applied to the mock catalogue (i.e. no 'missing galaxy' 
has been added to fill in the gaps), as well as all the other 
characteristics of the survey (e.g. the weight associated to 
galaxy shape measurements is also preserved). 

(iii) Following the procedure described in Section 2, the 
mass reconstruction is performed on a regular 512 x 512 grid 
which is a 2 X 2 re-binning of the simulations native pixel 
grid. The mass reconstruction is performed for 5 different 
smoothing scales, from approximately 2 to 9 arcminutes in 
radius. 

(iv) For some of the following tests performed on sim- 
ulations, we will be using the noise-free convergence map, 
obtained from the stacking of convergence maps at differ- 
ence source redshifts using the redshift distribution from 
(ii) above. These noise-free maps are maps directly obtained 
from the simulated Hght cone, without any connection to the 
mock catalogue. For this reason, and to make a clear distinc- 
tion with the quantity , we call Ksim the stacked noise free 
convergence map. 

Figure [l] shows a mass reconstruction example of one 
the tiles. The background image with the white contours 




Figure 2. For a smoothing scale of 6o = 2.5 arcmin, the different 
set of points show the cross-correlation profile {Ca-bli between 
two convergence maps Ka and rej , where a and b are one of " obs" , 
"B", "sim" or "ran". The black filled circles show {Coba;sim)r, the 
red filled triangles show {Cobs;s)r- Error bars for the filled circles 
and triangles are the Icr rms of the average over the 10 lines-of- 
sight. The light blue filled area shows the Icr region of (Cobs ran)r 
averaged over 100 random noise realizations. The Itr region in the 
latter also represents the deviation of the average. 



shows the reconstructed mass; the white contours represent 
the 1,2,3,4(7 levels, where a is the convergence rms deter- 
mined from the reconstructed map. The black contours rep- 
resent the 1, 2, 3, 4(T levels on the noise free convergence map 
Ksini- This figure illustrates qualitatively how realistic shot 
noise (ellipticity noise) noticeably affects the position, am- 
plitude, and sometimes even the presence of peaks in the 
reconstructed map. Many reconstructed peaks do not match 
a real mass peak and the opposite is also true. Although a 
quantitative analysis of peaks is left for another paper, this 
illustration is consistent with previous peak statistic anal- 



ysis such as Van Waerbeke (20001; Athreya et al. (20021, 



which shows that individual peaks are relatively noisy ob- 
jects, with a higher chance of being a coincidence (20 per 
cent for a 3(j peak) than Gaussian statistics. Masks are also 
shown in Figure [T] as the black areas with sharp boundaries. 
It shows that our mass reconstruction procedure does not 
generate catastrophic edge effects near the masks. The rest 
of this Section explores the reliability of the mass recon- 
struction quantitatively. 



3.3 Analysis of the mock catalogues 

In this Section we use our simulated data to verify that the 
level of shape noise and masks in the CFHTLenS data will 
not introduce systematics errors in the statistical properties 
of the reconstructed convergence map, such that the cosmo- 
logical information is left intact by the mass reconstruction 
procedure. This will be performed through two distinct tests. 

For the first test we focus on the cross-correlation map 
between two maps Ha and k;,, where a and b can be any of the 
signal Kobs , noise- free Ksim , galaxy rotated k± or randomized 
Kran maps. The cross-correlation map Ca\b is given by: 



Cat 



(17) 



where {...)o denotes the zero lag value of the auto-correlation 
map. This definition guarantees that Ca;f, is normalized, i.e. 
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the central pixel ofCa;b is equal to one if Ha = Kb- Ca-b is then 
azimuthally averaged within annuli as function of distance 
r from the central pixel; this quantity is called (Ca;6)r- 

Figure [2] shows {Ca;6)r for different combinations of re- 
constructed and noise-free maps for the Gaussian smooth- 
ing scale 00 = 2.5 arcmin. The filled triangles show {Ce:±}i, 
the cross-correlation profile between the reconstructed mass 
map and the galaxy rotated map, averaged over the 10 lines- 
of-sight of the mock catalogue (see Section 3.2). The error 
bars represent the error on the average which show that 
(Cobs;±)r is conslsteut with zero. The solid area shows the 
scatter of {Cobs;ran)r average over 100 random realizations 
of Kran. This quantity is also consistent with zero, which 
confirms that reconstructing the mass from random galaxy 
orientations does not generate systematics associated with 
for example edge effects at the field boundary or masks. The 
filled circles show the cross-correlation profile between the 
mass reconstruction and the noise-free map {Cobs;sim)r. The 
error bars illustrate the dispersion of the average over the 10 
lines-of-sight; the zero-lag cross-correlation coefficient is only 
0.3 which is another illustration that mass maps are noisy, 
especially for small smoothing scales. Larger smoothing 
scales show stronger cross-correlation coefficients, but also a 
larger correlation length. This first test demonstrates the ab- 
sence of systematic peaks around masks or low-density areas, 
if these were present we would expect a significant non-zero 
residual cross-correlation between the reconstructed mass 
map Kobs and k± or the pure noise reconstruction. 

The second test consists in looking at various higher- 
order moments of the convergence. Similar to the traditional 
measurement of the second order moment of the shear {7^) 
(here =71+72), we are interested in the true moments 
of the convergence {k"^), where n = 2,3,4,5. A measure- 
ment of the convergence moments from noisy reconstructed 
mass maps has never previously been reported. There is a 
complication with the measurement of convergence moments 
however: the convergence can only be measured within a 
smoothing window, as opposed to the shear which is given 
for each individual galaxy. Therefore the convergence noise 
is correlated for different points on the grid. In other words, 
the moments of Kobs and ke are not equal. A shot noise re- 
moval or de-noising procedure is therefore required in order 
to remove this bias from the convergence moments. In ad- 
dition, residual systematics tests can be performed on the 
rotated galaxies 'mass map' kb, and the reliability of mass 
maps can be evaluated by comparing the moments of to 
the moments of the noise- free maps Ksim ■ 

We choose a direct approach for the de-noising proce- 
dure, which consists of removing the contribution caused 
by the correlated noise from the measured moments. Fig- 
ure [2] shows that there is no cross-correlation between Kobs 
and the noise Kran; the absence of correlation means that, 
to first approximation, the signal and noise can be treated 
as statistically independent, as discussed in [Van Waerbeke] 
( |2000[ ). This lack of correlation can be used to derive rela- 
tively straightforward relations between the moments of Kobs 
and Ke- For instance, the observed second order moment 
(Kobs) is the quadratic sum of {k%) and (K?an)- Assuming we 
have a large number of random noise mass reconstructions 
Kran, the average moments of the noise map can be measured 
and subtracted from the observed signal moments. The re- 
lations below show how, for each smoothing scale 60, the 



observed moments of Kobs are related to the true moments 
of K_B and the moments of the noise map Kran: 



(KE)eo — («^obs)g^ — (K?an)eQ 

{ke}0o = («^obs)g^ 

{Ki)0o = {Kths)g^^ - 6(«:obs'^?an)9^ - (kLu) 

(4)eo = ('^obs),^ - 10(^L'^?an),„. (18) 

For clarity, (■••)fl(, denotes the moment average over one map, 
{---)0g is the same as (•••)eo ^^t averaged over several noise 
maps Kran, and {---)eo is the de-noised moment. We limit 
our analysis to higher order moments up to the fifth order. 
Beyond this we find the measurement to become very noisy 
with the CFHTLenS data. The Appendix includes the ex- 
pressions of the residual systematics moments, i.e. all the 
combinations of E and B modes which we expect to vanish 
if residual systematics are negligible, for n = 2, 3, 4, 5. As the 
measurements will show, the residual systematics are indeed 
consistent with zero. Our moment de-noising approach auto- 
matically takes care of varying noise across the field because 
noise variations for different positions on the maps are pre- 
served for the different noise realizations. A further improve- 
ment of the de-noising technique was implemented by tak- 
ing into account the fraction of pixels which are not masked 
within each smoothing window. This is achieved by con- 
structing a filling factor map computed from the smoothed 
grid with pixels initially set to one and then masked and 
smoothed with the same window applied to the lensing data 
for the mass reconstruction. In this way, regions around 
masks are down-weighted without biasing the signal. In or- 
der to eliminate the noisiest regions, a cut-off of 50 per cent 
is applied, and pixels below this threshold are excluded from 
the analysis. We verified that the application of this cut-off 
has a negligible effect on the measurement with a cut-off 
varying from zero to 80 per cent. Of course, large cut-off 
values boost the noise and it becomes difficult to evaluate 
changes in the signal. This procedure guarantees that pixels 
in the mass map contribute proportionally to the unmasked 
area located in a smoothing window. In the following, the 
de-noised moments of Kobs will be compared to the noise free 
moments measured on Ksim. 

Figure [3] shows the convergence moments as measured 
from the mock catalogues using the de-noising technique de- 
scribed here. We see that the de-noised moments (k^) (filled 
circles) are consistent with the noise-free moments (k";^). 
Note that the noise-free moments are perfect in the sense 
that they are directly computed from the Nbody simula- 
tions, there is no mass reconstruction, therefore no effect of 
the noise or discrete sampling and there is no mask. On the 
other hand, the moments (kJ) are de-noised moments ob- 
tained from the noisy mass reconstruction maps that include 
masks. The two moments agree which demonstrates that 
our mass reconstruction procedure leads to reliable mass 
maps where the statistical and cosmological information is 
preserved. On Figure [S] the de-noised moments errors were 
computed from the variance between the 10 lines-of-sight 
and divided by y/W. The noise-free moments error is shown 
by the filled region around the solid line, and it is also di- 
vided by vTO. The open symbols on Figure [3] show, for the 
different moments order, the residual systematics obtained 
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2, 3, 4, 5, measured on the simulations. Filled circles show the de-noised moments for 
the reconstructed mass map Kobs- The error bars show the \a deviation of the mean obtained from the rms of the 10 line-of-sight. The 
solid line inside the light blue region shows the noise-free moments measured on n^x^, and the filled light blue area shows the Icr deviation 
of the mean of the average 10 noise- free maps. Open symbols show all the possible de-noised combinations of Kobs and kj^: (K^)eo for 
the top left-panel, (K^)gg, (KEK%)eo' {'^%i^B)eo for the top-right panel, {k'^)bo, {i^E'^%)eo, {'^%'^%)eo' ('«|;'«s)9o ^^'^ bottom-left 
panel and (k;^)^^, {KEK,%)eo, (k|;K^)so, {i^%i^%)eo, {i^%f^B)eo for the bottom-right panel. 



from all the possible combinations of the reconstructed maps 
Kobs and K±, where the de-noised moments were computed 
using the expressions derived in the Appendix. With no ex- 
ception, all moments involving one or more of the rotated 
galaxies maps kx are consistent with zero, showing the the 
_B-mode is consistent with pure noise. 

These results on mock catalogues validate our approach 
for the CFHTLenS data and demonstrate that for the 154 
deg^ of the survey our reconstruction process is stable. Fur- 
thermore we have shown that realistic masking geometry 
does not alter the reconstruction, and therefore we are able 
to reliably quantify some of the most basic statistics of the 
projected mass density. 



3.4 Beyond the moments: convergence PDF and 
peak statistics 

The first few convergence moments give only a partial de- 
scription of the histogram of the convergence, also called the 
1-point Probability Distribution Function (PDF). However, 
the 1-point PDF, along with its extensions to higher order, 
contains additional information about the moment hierar- 
chy, and therefore about the structure formation process. 
It has been demonstrated for instance that a combination 
of different moments of the PDF, through the Edgeworth 
expansion, probes different aspects of the gravitational col- 



lapse as described in'Juszkiowicz et al.VlQQSV General char- 
acteristics of the PDF such as the height and the min- 
imum K of the convergence histogram are two important 
features which are not easily captured by convergence mo- 
ments. Therefore, it would be interesting to use the conver- 
gence PDF itself as a cosmological probe. This possibility 



has been theoretically explored (Munshi et al. 20041 using 



the aperture filtered shear, the same filter which transforms 
the shear 7 into a local, scalar quantity. Doing this analy- 
sis on the convergence k PDF would exploit the fact that 
the long wavelength modes are preserved, and presumably 
boost the signal-to-noise; such an analysis, based on the 
CFHTLenS mass maps, is left for a forthcoming study. 

Nevertheless we can already illustrate with our simula- 
tions the expected level of 1-point PDF signal-to-noise. Fig- 
ure|4]shows the measured 1-point PDF for k compared to the 
Kx 1-point PDF. The average 1-point PDF obtained from 
pure noise reconstructions (with error bars) is also shown, 
and the B-mode PDF is consistent with pure noise. This 
result is in agreement with Figure [S] that shows a negligible 
-B-mode for the convergence moments. One can also clearly 
see that the cosmological signal broadens the PDF compared 
to pure noise (or B-mode) PDF's. This 1-point PDF study 
opens up a new statistical analysis opportunity for gravi- 
tational lensing surveys which remains to be exploited. Of 
interest is the construction of alternative probes for cosmol- 
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Kran where the galaxy orientations have been randomised before 
the mass reconstructed. 



ogy using mass maps. For instance the minimum value of the 
convergence found in voids could be used as a direct mea- 
surement of the mass density parameter f2o that is relatively 
insensitive to the non-linear clustering. 

Another interesting use of convergence maps is the 
study of peak statistics. The idea of using peaks in dark 
matter maps to constrain the halo mass function was first 
proposed in [Van Waerbeke| ( |2000[ ) ; |Jain fc Van Waerbeke| 
(20001, it strongly relies on the properties of peak statis- 



tics in gaussian noise developed for the cosmic microwave 
background ( [Bond fc Efstathiou||1987[). The peak sta tistic 
has been recently pursued further ( Maturi et al.|20To Mar- 



ian et al. 



20121 



Peak statistics provide a higher level of 
statistical analysis of maps, where the sensitivity to cosmol- 
ogy and residual systematics is expected to be different in 
comparison to the more traditional moments statistics. It is 
important to distinguish between the maps constructed in 
this paper and the maps being discussed in [Marian et al.| 
( 2012[) which use the aperture mass filter, originally defined 
in |Schneider et ah] ( |1998[ ) . The aperture mass filter is a pass- 
band filter as opposed to the Gaussian (or top-hat) smooth- 
ing window used in our study, which is a low-pass filter. The 
low-pass filter preserves the large scale modes, hence keeping 
visible the voids and large overdensities the pass-band filter 
(see Section 4.4 below). The peak statistics of a low-pass 
map allows the study of the dark matter halo mass func- 
tion as function of the larger scale dark matter environment, 
which is essentially not possible with a band-pass filter. The 
two approaches are complementary and would benefit in the 
future being unified in a single approach of mass map mak- 
ing. Peak statistics, convergence PDF and the morphological 
analysis of large scale structure from the CFHTLenS data 
are left for future studies. 



4 RESULTS FROM CFHTLENS 

Residual systematics in the shear signal, based on the se- 
lection of "good" fields passing the systematics tests, are 
given in|Heymans et al.| ([20T2b|). In total, 129 MEGACAM 



pointings (out of 171) passed the residual systematics tests. 
The convergence moments presented in this work are based 
on the good fields list, while the mass reconstruction is per- 
formed on all the fields. Each galaxy in the CFHTLenS cat- 
alogue has a shear estimate e°''^ and a weight w tMiUer] 
et al.|20T3 l, a calibration factor m ( Heymans et al.||2012b l, 
and a photometric-redshift probability distribution function 
( Hildebrandt et al.||2012 l. In the following Sections we de- 
scribe how these elements are integrated into the mapmak- 
ing process. 

For the map making, we select galaxies in the redshift 
range z = [0.4, 1.1], which is sufficiently broad to avoid the 



complication of intrinsic alignment (Heymans et al. 20131, 
but narrow enough to guarantee that the 3-point signal is 
detectable at a few sigma ( |Vafaei et al.|2010| ). 



4.1 Mass reconstruction 

Mass reconstruction on the CFHTLenS data follows the 
same procedure as described in Section 2 using the KS93 



algorithm (Kaiser & Squires 19931. The regular grid over 



which the mass reconstruction and smoothing are performed 
consists of square pixels with an area of approximately 1 
arcmin^. The average of the galaxy ellipticities is first cal- 
culated within each pixel, using the calibration correction 



implemented in Miller et al. (20131 in order to account for 



the constant correction C2 (Heymans et al. 2012b I, the multi- 



plicative shear measurement bias (l-t-rrii) and the weighting 
Wi for each galaxy located at Qi within that pixel. The shear 
estimate per galaxy is not divided by (1 -I- m^), but instead 
the average epix in each pixel is given by the weighted sum 
over the galaxies in that pixel: 



(19) 



The smoothing of this grid-averaged ellipticity map is then 
performed using the Gaussian window function given by 
Eq.([9|. Mass maps are reconstructed with the following 
smoothing scales 1.8, 3.5, 5.3, 7.1 and 8.9 arcmin. The fi- 
nal maps have a size of 512 x 449, 512 x 505, 512 x 495 and 
512 X 502 for Wl, W2, W3 and W4 respectively, which is 
chosen to be close to the grid size that was used on the sim- 
ulations. The fact that the angular resolution is different for 
the different fields has no effect on our results. Mass recon- 
struction is also performed after rotating galaxy ellipticities 
by 45 degrees in order to probe _B-modes. The analysis of 
the cross-correlation map following Eq.( |17| l shows that the 
B-mode in the data is consistent with zero (not shown here). 

4.2 Cosmic statistics on CFHTLenS mass maps 

Following the procedure described in Section 3.2, 100 in- 
dependent pure noise maps Kran are reconstructed for each 
CFHTLenS field and each smoothing scale. The pure noise 
maps are necessary for the de-noising of the convergence 
moments. This paper is the first to measure cosmic shear 
statistics directly from reconstructed mass maps and it is 
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Figure 5. Redshift distribution of the galaxies from the 
CFHTLenS catalogue that are used for the measurement of the 
convergence moments. Data points show the actual probability 
distribution obtained from the photometric redshifts and the solid 
line shows the best-fitting function defined in Eq |20| The best- 
fitting parameter values are a = 1.50, b = 0.32, c = 0.20 and 
d = 0.46. 



therefore it is important to quantify the reliability of the 
measurements against predictions. Theoretical predictions 
for the second and third order moments are calculated us- 
ing Eq ^l4| and |17[ The key ingredient for these predictions 
is an accurate redshift distribution 713(2) of the galaxies. 
To estimate ns{z) we use the probability distribution func- 



tion of the photometric redshifts (Hildebrandt et al. 20121 
whose robustness has been thoroughly tested in |Benjamin| 
et al. (20121. The implementation of this approach in other 



CFHTLenS papers, where the second order moment of the 
shear has been used to constrain cosmological parameters, 
shows that it tends to a consistent interpretation of the weak 



lensing data as a function of source redshift ( Benjamin et al. 
20T2i IHeymans et al.|[20T2a| [Kilbinger et al.||2013| |Simpson 



et al. 12013 1. 



The data points in Figure |5] show the stacked redshift 
distribution PDF for our galaxy selection with photometric 
redshifts 0.4 < z < 1.1. In order to make predictions for the 
lensing statistics, we fit the observed redshift distribution 
with the following four-parameter double- Gaussian model: 



ns{z) — axexp 



b2 



-l-cxexp 



d2 



.(20) 



Figure |5] shows that this is a reasonable fit to the data, in 
particular in capturing the long, slowly-decreasing, tail at 
high redshifts. The best fit values are a = 1.50, 6 = 0.32, 
c = 0.20 and d = 0.46. 

Figure [6] shows the second to fifth order moments of 
the convergence. The solid lines show the predictions for the 
second and third order moments using the redshift distribu- 
tion derived above. The predictions are not a fit to the data. 
They are computed using the CAME non-linear power spec- 
trun0for the ACDM model = 0.0471, ^cdm = 0.2409, 
Qa = 0.712 for the baryonic, dark matter, and cosmological 
constant density parameters respectively. The matter power 
spectrum is normalised to erg = 0.792 at 2 = 0. The other 



http:/ /camb.info/ 



cosmological parameters such as the power spectrum tilt, 
running spectral index, etc. are identical to the best-fitting 



values from WMAP7 (Komatsu et al. 20111. This partio 



ular choice of parameters is consistent with the parameter 



values obtained in the CFHTLenS companion papers ( 


Ben- 


jamin et al.||2012| |Heymans et al.||2012a| 


Kilbinger et al. 


2013| Simpson et al.| 2013| Kitching et al. 


20131). The 


pre- 



dictions agree well with measurements for the second and 
third order statistics, which gives strong confidence in the 
reliability of mass maps to extract cosmological information 
and it confirms the analysis performed on the shear. The 
error bars represent the sampling variance between the four 
CFHTLenS fields, accounting for the different image sizes, 
therefore they capture the total error budget (including sta- 
tistical variance). There are no reliable fully non-linear pre- 
diction for the fourth and fifth order moments, therefore a 
full cosmological analysis including all moments must rely 
on ray-tracing simulations. This work is left for a future 
study where the properties of the full convergence proba- 
bility distribution function will be investigated. As can be 
seen in Figure |6j the CFHTLenS data show a marginal de- 
tection of the fourth order moment and no detection of the 
fifth. The residual systematics are consistent with zero for 
all moments. 



4.3 The connection between Large-scale dark 
matter and baryons 

Sections 3 and 4 describe the mass reconstruction process in 
the presence of realistic noise and masks. A cosmological sig- 
nal is clearly measured on the CFHTLenS data and the level 
of residual systematics is shown to be consistent with zero. 
In this section we are interested in the comparison between 
the reconstructed dark matter and the matter distribution 
of the stellar content. The connection between dark matter 
and baryons is quantified in the CFHTLenS companion pa- 
pers ( [Gillis et aL]|2013| [Velander et al.||2013[ ) at the galaxy 
and galaxy-group scales. Here we are interested in the con- 
nection between dark matter and baryons at much larger 
scales, exploring if there are large scale features that extend 
out to, for example, a void size which is common to the dark 
matter and the stellar content. To this end, we develop a new 
approach which first consists in predicting the dark matter 
map from the stellar mass distribution ( Wilson et al.|200T |, 
and then compare the peak distribution from the predicted 
K map to the lensing map. 

In order to construct a predicted k map, a Navarro- 



Frenk- White dark matter density profile (Navarro et al. 
|1997[ ) is assigned to each galaxy in the CFHTLenS cata- 
logue; the procedure is detailed below. There is no need to 
separate the galaxy population into sources and lenses be- 
cause most galaxies act as both due to our broad redshift 
selection. Galaxies are assigned a dark matter mass follow- 
ing the relation between stellar mass and halo mass from 
[Leauthaud et alT] ( |2012[ ) , which provides a relation for red- 
shifts between 0.2 and 1, covering almost the entire redshift 
range of galaxies in our work. The only complication is that 
this relation is only provided for central galaxies, while in 
our sample a large fraction are satellites. Unfortunately it 
is impossible to separate the central galaxies from the satel- 
lites on an individual basis, it can only be done statistically 



(Velander et al. 20131. We therefore proceed by assuming 
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Figure 6. Moments of the convergence (k"), n = 2,3,4, 5 measured on the CFHTLenS data. Error bars show the la deviation from the 
mean of the the four CFHTLenS fields. Solid lines are the moments measured from the signal maps and de-noised using the procedure 
described in Section 3.2. Open symbols show the different de-noised combination of the signal map K^hs a^nd systematics map similar 
to that shown in Figurejslfor the simulations. The solid line shows the second order moment (top-left) and third order moment (top-right) 
predictions from Eq^ 14| and | f 7| using the WMAP7 cosmology (see text in Section 4.3). 



that all galaxies are central galaxies and all follow the stel- 



lar mass to halo mass relation in Leauthaud et al. (20121. 



This will overestimate our predicted total mass, but to a first 
approximation it should not dramatically affect the relative 
distribution of mass. Using the Figure 4 from |Leauthaud| 
et al. ( 2011 1 we anticipate that, on average, there are roughly 



one to two satellite galaxies for every central galaxy. This 
should lead us to overestimate the total mass by roughly a 
factor of 2 to 3. The exact calculation is not needed as we 
are only interested in an order of magnitude estimate of how 
wrong our predicted convergence can be. In a future work, 
the same strategy will be applied to clusters instead of indi- 
vidual galaxies, which should mitigate this effect. In order 
to complete our convergence prediction from the galaxies, 
we need to assign a concentration to each halo. To this end, 
we use the mass-concentration relation calibrated from nu- 
merical simulations in |Munoz-Cuartas et al.| ( |2011[ ) . 

At this stage, each galaxy in the CFHTLenS catalogue 
is associated with a dark matter halo of known concentra- 
tion and mass. The last step is to apply the lensing kernel 
in order to predict the convergence based on the redshifts 
of the lenses and sources. Every galaxy is simultaneously 
both a lens and a source, depending whether they are in the 
background or the foreground relative to other galaxies. One 
can then compute for each galaxy a predicted convergence 
based on the foreground mass distribution coming from all 
galaxies located at lower redshift. For a source galaxy at 



location 6 on the sky with redshift zs, and A'' foreground 
lenses at redshifts zl^ the total convergence predicted from 
the baryonic distribution is given by: 



Kgal(0) = 



JV 



E,(|6>-6/.l 



it(ZL,,Zs) 



(21) 



where T.i{\6 — Oi\) is the projected halo mass of lens i cen- 
tred at 9i, and T,crit{zLi, zs) is the critical density given 
by: 



TicritizLi , Zs) 



fK{ws) 



AnG fK{wL)fK{ws — WLi)' 



(22) 



Note that the average predicted convergence Kgai is cal- 
culated only after all haloes have been assigned to the 
galaxies. The critical density depends on the observer-lens, 
lens-source and observer-source angular diameter distances 
/k(wz,), fK{ws - WLi) and /x(w^s). The sky-average pre- 
dicted convergence is set to zero by subtracting the mean 



Kgai in Eq. (21|. We assume that the weak lensing approxi- 



mation applies which means that the convergence from the 
different lenses can be added linearly. It is important to em- 
phasize that the lens redshift goes down to the lowest reli- 
able value z = 0.2, and that the sources only cover redshifts 
z = 0.4 and higher (in order to be consistent with our source 
galaxy selection outlined in Section 4.2). 

The convergence predicted from the baryonic content 
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Figure 7. In the Wl field, comparison of the total mass from gravitational lensing reconstruction (underlying continuous map with 
contours) and the predicted mass inferred from the foreground galaxies is indicated by the open white circles in the map. The field-of-view 
is approximately 9x8 deg^. Open circles indicate the position of peaks with positive curvature in the predicted projected mass. The 
number density of circles is therefore indicative of the amount of mass expected in projection in a particular region. The size of the 
circles are also indicative of the height of the peak of the predicted mass map is. The open circles distribution coincides with the high 
density regions in the reconstructed mass density map. 



is assigned to each galaxy in the CFHTLenS cata- 
logue. Following Section 4.2, the Kgai(S) are placed on the 
same regular grid that is used for the lensing mass recon- 
struction. The lens&t weighting is used to determine the 
average «:gai(6') within each pixel, and the same Gaussian 
smoothing is applied. For most galaxies, the size r2oo of a 
single galactic halo is comparable to the size of a pixel, which 
is of the order of half an arcminute, in the final resolution 
map. The galaxies in the masked regions therefore do not im- 
pact on our analysis as the extension of their halo is small 



The Gaussian smoothing window is sometimes truncated by 
the masks, but this effect is minimal, because of the filling 
factor cut of 50 per cent applied to the predicted mass map. 
Furthermore, this is a random, zero net effect, leading only 
to larger noise around the image masks, but not to a bias in 
the projected convergence. 

In order to compare the reconstructed lensing map and 
the predicted map, we adopt the following strategy: positive 
curvature peaks are first detected in the Kgai map, with a 



Note that the haloes from low rcdshift galaxies do extend over 



several pixels, but their lensing 
proximity to the observer 



efficiency is small due to their 
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Figure 8. Left panel: The x-axis values show the predicted convergence k obtained from the galaxy distribution taken at positive 
curvature peak locations for all W fields. The conversion of baryonic mass to total matter convergence is detailed in Section 4.4. The 
y-axis shows the convergence k measured from the mass reconstruction map using gravitational lensing. Individual small dots on the 
figure represent a peak in the predicted k map, and the lensing k is taken from the same location. The black open circles show the binned 
average, and the red open circles show the binned average on the 45 degrees rotated galaxies _B-mode reconstructed map. Error bars 
represent the dispersion on the mean. Right panel: now the x-axis (predicted k from galaxies) is taken at negative curvature peaks. The 
detection of voids is manifest in the very significant negative values of the lensing reconstructed k for negative values of the predicted 
convergence. 



peak being defined as a pixel where all the pixels around 
have a lower amplitude. A signal-to-noise (or significance) is 
then associated to each peak with the noise being k r.m.s. 
in the lensing map computed at the peak location using the 
pure noise reconstruction. Figure [7] shows the location of 
Kgai peaks for Wl (shown as white circles) against the recon- 
structed lensing map in the background. Contours are shown 
for the background at 1,2,3 and 4 sigma, which is also an 
indicator of the significance of structures in the lensing map. 
Note that the lensing map is shown for the largest smooth- 
ing scale of 8.9 arcmin, while the predicted peaks have been 
identified on the 1.8 arcmin smoothing map. On average, 
the number density of ftgai peaks coincides with the lens- 
ing mass overdensities. A quantitative comparison between 
the predicted convergence and the lensing convergence is 
shown in the left panel of Figure |8] the small dots in this 
figure show, for each peak with positive curvature detected 
on the Kgai map, the corresponding value of the lensing map 
Kiensing at the Same location, for the same smoothing scale 
of 1.8 arcmin. Although only positive curvature peaks have 
been used in this plot, it does not mean that the lensing 
convergence is positive. In fact, one can have haloes inside 
low density regions with negative convergence. One can see 
that this is the case, although on average positive curvature 
peaks live in positive convergence regions. A comparison of 
the binned Kgai peak values to the averaged lensing conver- 
gence in the bin shows that the reconstructed and predicted 
convergence are strongly correlated, consistent with the fact 
that baryons trace dark matter to first approximation. The 
lensing convergence is roughly 2 — 3 times smaller than the 
predicted convergence; this is expected, as discussed earlier 
in Section 4.4 because we assumed that all galaxies are cen- 
tral galaxies of the host dark matter halo. Following the 
numerical calculations from |Leauthaud et al.|p011[ ) this er- 
roneous assumption would lead to an overestimate of the 
predicted mass by a factor 2 — 3. Note that for an individ- 



ual peak, the noise is large with a r.m.s. of cr^ — 0.015, 
consistent with the fact that within each smoothing window 
there are approximately ~ 100 lensed galaxies. Figure Hjalso 
shows the peak analysis performed with the B-mode lensing 
map, and the correlation between kb and Kgai completely 
vanishes. 

The analysis of the positive curvature peaks provides 
a consistent picture in which the galaxy distribution is 
strongly correlated with the lensing mass, with more peaks 
and more mass in higher density regions. A similar anal- 
ysis can be performed on negative curvature peaks. This 
is statistically equivalent to probing underdense regions, or 
voids. Figure [9] shows such an analysis for Wl where white 
triangles represent negative curvature peaks (local minima) 
in Kgai- One can see that the triangles preferentially popu- 
late underdense regions of the lensing map; this is a con- 
vincing illustration of the detection of giant projected voids 
in large-scale mass maps based on lensing data. Compari- 
son of Figure [7| and [9] reveals how the positive and nega- 
tive curvature peaks probe large-scale structures and voids, 
respectively. Right panel on Figure [8] shows how negative 
curvature peaks in the predicted map Kgai compare to the 
lensing mass reconstruction Kicnaing. It shows that the de- 
tection of voids is significant, with residual systematics also 
being consistent with zero in the underdense regions. Nega- 
tive curvature peaks tend to be located in majority in nega- 
tive convergence regions, as one would expect. The predicted 
convergence Kgai is closer to the reconstructed convergence 
Kicnsing for the voids (Figure ??) and for positive curvature 
peaks (Figure[8|. A possible explanation would be that there 
is less satellite galaxies in voids than in overdense regions. 

Figures [lO] and [iT] show the same analysis for the three 
other CFHTLenS fields W2, W3, and W4. The peak and 
void statistics shown in Figure |8] and ?? include all four 
fields. The key message from the comparison of the baryon 
and dark matter maps is that dark matter maps trace large 
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scale structures reliably, but evaluating the significance of 
most of the individual mass peaks is challenging for two 
reasons: first, the shot noise is high and only very few peaks 
are above the secure 4a threshold. Second, the maps clearly 
show that the large scale structures are a source of noise 



• Our approach consisted of removing the noise bias when 
measuring the statistics of the reconstructed mass. A com- 



comparable to the shot noise (Hoekstra 2001 Hoekstra et al 



20111 



5 CONCLUSION 

This paper is the first quantitative cosmological analysis of 
mass maps reconstructed from lensed galaxy shapes. We val- 
idated our approach first using N-body simulations and then 
applied it to the CFHTLenS data. We find that the tradi- 
tional cosmological interpretation of lensing data on the con- 
vergence is quantitatively possible, but also that the scalar 
nature of mass maps additionally opens a new window for 
lensing analysis where global features (such as the morphol- 
ogy of large scale structures and the peak statistics) could 
be used to constrain cosmological models. 

Using N-body simulations, we have shown that the 
reconstruction process is stable, and the input cosmologi- 
cal signal is recovered accurately, despite the presence of 
masks, the relatively high level of shape noise and the non- 
Poissonian distribution of the background sources. Mass re- 
construction was performed and tested with the traditional 



KS93 algorithm ( Kaiser fc Squires_1993J . The testing of the 
full non-linear mass reconstruction at the same level of pre- 
cision is left for a future study. 

The application to the CFHTLenS data shows for the 
first time that windowed statistics of the convergence can be 
measured on mass maps. We find excellent agreement be- 
tween the second and third order moments measured on the 
reconstructed mass maps and the predictions for a cosmo- 
logical model determined from a shear correlation function 



analysis of the same data ( Benjamin et al.||2012 Heymans 
et al.||2012a| IKilbinger et al.||2013| ISimpson et al.||2013l). 



An attempt at measuring higher order statistics shows a 
marginal detection of the fourth order moment. For all the 
moments of the convergence, the residual systematics are 
found to be consistent with zero. In Section 4.4, we compare 
the reconstructed convergence to the predicted convergence 
using galaxies as tracers of dark matter haloes, where a halo 
is assigned to each galaxy. We have shown that the predicted 
and reconstructed mass maps are strongly correlated to each 
other. The maps reveal the existence of large voids in the 
projected dark matter distribution which can span regions 
as large as 3 — 4 degrees on the sky. 

We consider this paper to be a feasibility study which 
strongly suggests that future precision-cosmology on mass 
maps is possible. Cosmology with mass maps enables studies 
that are currently not possible with the shear -shape- infor- 
mation only: non-local structure morphology, peak statistics 
with long wavelength modes included, cross-correlation with 
other cosmology maps, and statistics of convergence proba- 
bility distribution functions. However, before mass maps can 
really be turned into a useful cosmological probe, several is- 
sues will have to be addressed. 

• A correct non-linear reconstruction will need to be im- 
plemented in order to be able to account for the most mas- 
sive structures such as clusters of galaxies. 



pletely different strategy is to de-noise the map itself ( Starck 
eraL][2006l ). It remains to be tested if the latter is a viable 
approach for precision measurements. 

• The photometric redshift uncertainty will need to be in- 
cluded in the dark matter-baryon comparison. How to con- 
struct a predicted mass map from a distribution of galaxies 
with photometric redshift errors is an open problem. 

• A potential fundamental limitation with maps is that 
the current mass reconstruction process cannot distinguish 
between intrinsic galaxy alignment, gravitational lensing 
or other systematic errors that lead to correlated galaxy 
shapes. For our particular choice of redshift range, and from 
the residual systematics studies, we know that our measured 
signal is mainly caused by gravitational lensing, but preci- 
sion cosmology will demand that we can identify and sep- 
arate the different causes. It is not currently clear whether 
this is possible with mass maps and further theoretical stud- 
ies will be needed. 
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Figure 9. Similar to Fig[7] witii the background map and contours mapping the total projected matter from gravitational lensing mass 
reconstruction, now the open triangles show the position of the negative curvature peaks in the mass map predicted from the foreground 
galaxies. The open triangles now unambiguously trace the underdense regions in the density mass map. Larger triangles correspond to 
deeper negative curvature peaks. 
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the science analysis and interpretation in this paper, the founding 
core team and those whose long-term significant effort produced 
the final CFHTLenS data product. The second group covers mem- 
bers of the CFHTLenS team who made a significant contribution 
to either the project, this paper, or both. The CFHTLenS collab- 
oration was co-led by CH and LVW. 



6 APPENDIX 

For the measurement of moments from maps, it is necessary 
to remove the contribution from the noise component. It 
is not possible to remove the noise at the pixel level on 
the map, but it is possible to remove it statistically from 
the noisy moments. Our approach is to assume that, for a 
smoothing scale 6o, the observed convergence map Kobs is 
the sum of the true signal ke and an uncorrelated noise 
component Kran- We are repeating Eq ^lO| and|ll| 



(23) 



The observed B-mode convergence map k,±, which is the 
reconstructed mass from the 45 degree rotated galaxies, is 
similarly related to the true B-mode k.b and the noise con- 
tribution Kran: 



— Kb + 



(24) 



Note that the true B-mode should be zero if there is no 
residual S-type systematics. Since we are precisely testing 



this hypothesis, we keep as a quantity to be measured 
from the data. The noise contribution to each convergence 
moment measured in this work are then given by: 



(kek|)£)o = (Kobs^i)^^ 



I 2 2 \ _ i 2 2\ 7 2 2 ^ 7 2 2 V 

(,KeKB/9o — (^^obs^^i ^ V'^^obs'^ranjg^ — V'^±^ranjg^ 

— (Kran)9jj 



/ 2 3 \ / 2 3 \ 

3 2 \ 



9o 



(25) 



The left-hand side is the de-noised (unbiased) moment and 
the right-hand side shows the combination of moments lead- 
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ing to this unbiased measurement. The measured signal, B- 
mode, and random noise maps are respectively given by ^obs, 
K±, and Kian- As a reminder from the definition given in 
Section 4.2, a few ensemble averages have been introduced: 
{■■■)gg denotes the moment average over one map, {■■■)g^ is 
the same as {■■■)gg but averaged over several noise maps Kran, 
and {■■■)eQ is the de-noised moment. One can see from these 
equations that there are as many systematics moments as 
there are permutations of k_e and kb for a given moment 
order. A total of hundred random maps have been gener- 
ated for each smoothing scale. In order to derive Eqs.(25l 



we assumed that all the terms with odd order moments of 
the noise maps Kran can be neglected as they only introduce 
noise in the equation above. 
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